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We discuss the effects of a novel polaronic disorder in the recently proposed two-fluid model 
for manganites. Using effective field theory as well as direct numerical simulations, we show that 
this disorder can have dramatic effects in terms of the transition from ferromagnetic insulator to 
ferromagnetic metal upon hole-doping, including an Anderson localized regime where variable range 
hopping may be observed. 



I. INTRODUCTION 

Doped rare-earth manganites Ai-^AE^BOs (A=Rare- 
Earth, AE=Alkaline Earth, B=Mn) are transition metal 
oxides where there is an intricate interplay of charge, or- 
bital, spin and phonon degrees of freedom. No less impor- 
tant, however, is the role of disorder, which has dramatic 
effects upon the charge ordering, magnetic transition, 
etoi. While many different types of disorder has been dis- 
cussed in the context of manganites (eg. A-site disorder, 
B-site disorder etc.)2~— , in this paper we will introduce a 
novel kind of disorder apparently present in manganites 
but not widely talked about, and demonstrate its effects 
using the two-fluid ('£-6') model Hamiltonian 5 which was 
proposed earlier for manganites. 

The basic physics of manganites involve a coexistence 
of fast moving band electrons, and self-trapped JT po- 
larons, interacting by Coulomb interaction, and strongly 
coupled to a background of core spins. The l l-V model 
is an effective low energy Hamiltonian which implicitly 
captures the crucial effects of these interactions and the 
quantum dynamics of the JT phonons. It invokes two 
types of e g electrons, one polaronic and localized (£), and 
the other band-like and mobile (6), and is given by, 
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The polaronically trapped l t species has site energy 
—Ejt (~ 0.5 eV), and an exponentially reduced hop- 
ping (~ 1 meV) which has been neglected, while the 
non-polaronic '6' species (site energy 0) has undimin- 
ished hopping t ~ 0.2 eV. Udd (~ 5 eV) is the ef- 
fective on-site Coulomb repulsion between the '£' and 



'6' species and Jh the aforementioned Hund's coupling. 
Jp is a novel ferromagnetic virtual double- exchange cou- 
pling (~ 10 meV) between the core spins, which arises 
from virtual, fast (adiabatic) hopping processes of the 
'£' electrons to neighboring sites and back, leaving the 
local lattice distortion unrelaxed^. The chemical poten- 
tial, fj,, imposes the doping determined filling constraint: 
E,(W + (n&o-)) = (1 - At T=0, in the fully po- 
larized ferromagnetic phase, it reduces to a form very 
similar to the Falicov-Kimball model 6 . 

While disorder is commonly grouped into 'annealed' 
and 'quenched', we make a further distinction in the con- 
text of this £ — b model for manganites, namely: intrinsic 
and extrinsic. We note that from the basic premises of 
the I — b Hamiltonian, the t electrons are polaronically 
trapped at specific lattice sites, and have no quantum 
dynamics. They however, equilibriate according to clas- 
sical statistical mechanics. The mobile b electrons hop 
in the background provided by these static £ polarons. 
Hence, the £ species represents an annealed disordered 
background to the mobile b species. We refer to this dis- 
order as 'intrinsic', since it is 'self-generated' by the sys- 
tem under appropriate conditions. On the other hand, 
disorder due to doping at the A-site by alkaline earth 
atoms, or doping at the B-site by Aluminium, Chromium 
etc., is referred to as 'extrinsic'. Thus, the b electrons face 
two kinds of disorder of the compositional type: namely 
intrinsic, annealed disorder due to the I polarons^, and 
extrinsic, quenched disorder due to foreign dopant atoms. 
The effects of quenched disorder of the extrinsic type has 
already been considered in some detail in the literature, 
both at the abstract level^, and in the specific context 
of manganite o 10 ' 11 . In this paper we discuss the effects 
of this novel, intrinsic disorder in the context of the £ — b 
Hamiltonian for manganites. 

The £ — b model has been studied using Dynamical 
Mean Field Theory (DMFT) in Ref 5 , which is exact in 
the limit of large dimensions. Due to the specific sym- 
metries of the £ — b model, discussed in Ref£, it becomes 
exactly solvable in this limit, enabling a complete study 
of the phase diagram. However, DMFT, like its coun- 
terpart CPA (Coherent Potential Approximation), does 
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not include effects of Anderson localization which is a 
crucial effect of disorder at low temperatures, and which 
we want to focus on in this paper. The reason for this 
is discussed in more detail below. Hence, in this paper, 
we study the I — b model using an approach pioneered by 
Vollhardt and Economou et al which can take us beyond 
the DMFT approximation. We also compare our results 
with studies of Inverse Participation Ratio (IPR) done 
using a direct numerical simulation. 

II. DMFT OF L-B MODEL: SUMMARY 



where TL a p(uj,q) is the well-known current-current corre- 
lation function. In the DMFT limit of large dimensions, 
only the elementary particle-hole Drude bubble for this 
correlation function survives, while all higher terms in- 
volving the irreducible vertex become negligible — . The 
contribution to the optical conductivity for a hypercubic 
lattice then becomes^: 

1 1 d - 

a(iu n ) = — - 4sin 2 (ki)G(k, iv n )G(k, iv n +iuj n ) 

kv n o l — l 

(7) 



In the l—b Hamiltonian, the mobile b electrons move in 
an annealed disordered background of site-trapped I po- 
larons. The DMFT approximation can be used to study 
the electronic and transport properties of this Hamilto- 
nian. We first solve the DMFT exactly using semicircular 
DOS, and in the limit Udd~ > oo^, in the fully spin- 
polarised ferromagnetic phase at T=0i2. The DMFT 
approximation reduces the complicated many-site prob- 
lem given by Eqn [1] into an effective single-site problem, 
represented by the action: 

A D MFT{ni) = y~][(Gto?t h (iVn,toz) - Uddnijbnbn] 
n 

- {Ej T - n)nt (2) 

where np is the fractional occupancy of ^-electrons at the 
DMFT site, and G(oS) is the Weiss or bath Green's func- 
tion. Solving this single site 'impurity' problem, we ob- 
tain the full impurity Green's function G(uS). Since this 
action is similar to that for the Falicov Kimball model 6 , 
hence it can be solved exactly. We can also obtain the 
exact self-energy for the b-electrons using the Dyson's 
equation: 

E(w) =Q- 1 {uj)-G- 1 {lo) (3) 
This expression is given by: 

^)=(l (4) 

where 

g( u ) = - I (5) 

u + - W D 2 

and Wq is the probability of there being no £-polaron on 
the impurity site of the DMFT, to be determined self- 
consistently. For low dopings, Wq ps x. D is the half- 
bandwidth of the bare semicircular DOS. 

The optical conductivity of a thermodynamic system 
can be obtained within linear response, using the Kubo 
formula: 

1 d 

Ka(uj + ,q^ 0) = -— ^in aa (uj + ,q -> 0) (6) 

a — 1 



III. BEYOND DMFT: SELF-CONSISTENT 
THEORY OF LOCALIZATION 

The DMFT approximation does not capture the 
physics of Anderson Localization, which may well be im- 
portant for the low temperature transport behaviour of 
manganites having both intrinsic and extrinsic disorder. 
In particular, the Metal Insulator Transition (MIT) as 
a function of doping at low temperature is expected to 
be strongly influenced by localization physics. To inves- 
tigate this aspect, we go beyond DMFT and study the 
mobility edge behaviour as a function of doping using the 
Self- consistent Theory of Localization (STS) developed 
by Vollhardt et. al. 1A . Economou et. al. had proposed 
a particularly simple, and practically useful prescription 
for finding the mobility edge by mapping the STS equa- 
tions formally to equations for bound state formation in a 
potential well^. We use Economou's approach in this pa- 
per, which we refer to hereafter as Potential Well Analogy 
(PWA). We first find the single particle Green's function 
for the 6-electrons, as well as the self-energy, using the 
DMFT approximation at zero temperature. The dc con- 
ductivity of the system is also found within the DMFT 
approach. Then we plug these quantities into the PWA 
equations to obtain the mobility edge at any particular 
doping. 

In order to incorporate localization physics, one has to 
add correction terms to the Drude conductivity, obtained 
from 'maximally crossed' diagrams in the theory of weak 
localization^ 6 -: 

, , 2e 2 /■V'-i dk 

where do is the Drude conductivity, l e i is the elastic mean 
free path, while ij„ is the inelastic mean free path. Using 
the Einstein relation between conductivity and diffusiv- 

ity, 

a = 2e 2 pD (9) 

one can rewrite eq|8]purely in terms of diffusivites. Now, 
the theory was made self-consistent by Vollhardt and 
Wolfle 1 - by replacing Dq by the full D(ui) on the de- 
nominator of the integrand. 
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The localization length A was found by identifying 
tj^ — » ^j- in the insulating regime. 

Economou and Soukouli*^ pointed out that this equa- 
tion has the same formal structure as the equation which 
determines the decay length of a bound state in a poten- 
tial well. 
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where fl is the volume of the primitive unit cell, fcf, is the 
inverse decay length of the bound state, while \Vq\ is the 
well depth, and m* is the effective mass near the band 
bottom. 

For the case of tight binding on a hypercubic lattice, 
one particular site has the well potential Vq. Since the 
integral is over the Brillouin zone, one establishes a con- 
nection between the lattice constant a and the mean free 
path l e i, neglecting l/li n at low T. Now, knowing the 
minimum depth of the potential well necessary for the 
appearance of a bound stat o 7 ' 18 , one can estimate the 
mobility edge. 

The mean free path 1(E) is found as usual: 



1(E) = v(E)t(E) 

v (E - KT.(E))h 
2XT,(E) 



(12) 
(13) 



where Vq(E) is defined as the average speed of Bloch 
waves over the surface of constant energy E. 

The Drude conductivity (Jq(E) can be found either in 
its full glory from Eq [7] or using the weak scattering ap- 
proximation: 
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where So (E) is the area of the surface of constant energy. 



IV. RESULTS 



As discussed in detail in*, the I — b model shows a 
transition from Ferromagnetic Insulator to Ferromag- 
netic metal, as the doping is increased, just as observed 
in real manganites. This is because of depletion in the 
number of trapped Jahn- Teller polarons, and correpond- 
ing decrease in scattering and increase in paths avail- 
able for hopping, upon hole-doping. Within DMFT, for 
U — y oo, the critical doping for this Insulator-Metal tran- 
sition (IMT) can be obtained analytically as 
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which corresponds to the doping where the bottom edge 
of the &-band hits the sharp level at energy Ejt cor- 
reponding to Jahn- Teller trapped ^-polarons. Since for 
U — > oo the DMFT 6-electron DOS is proportional to the 
imaginary part of the bath Green's function mentioned 
above, they have the same bandwidth D = y/WaD, as is 
evident from Eq[SJ Hence, when the 6-band bottom hits 
the polaronic level, D = Ejt, giving the relation Eg [T5l 
above. For x < x c , the entire (1 — x) fraction of electrons 
were polaronically trapped in the local Ejt level, giving 
a polaronic insulator. For x > x c , the 6-band states are 
available for the electrons to delocalize, which according 
to DMFT, should immediately lead to a metallic state. 
We would like to improve upon this picture in this paper 
by bringing in the effect of Anderson localization. 

We calculate the relaxation time r, mean free path 
I and transport DOS v(E)S(E) for a cubic lattice and 
a single-orbital, tightbinding band. The DMFT parame- 
ters are chosen so that the bottom edge of the 6-band hits 
the sharp £-leve\ at a doping (x^ MFT ) of about 0.25. On 
account of the exact analytic relation mentioned above, 
this is accomplished by choosing, eg., Ejt to be —3 in 
units of hopping. This is because, for a cubic lattice, the 
half bandwidth is 6t, which we identify with the half- 
bandwidth D of the semicircular DOS mentioned above. 
Thus D = 6 in units of hopping t. The mobility edges 
are calculated both using the full conductivity, and the 
weak-scattering assumption as in Eq 1141 




x (doping) 

FIG. 1: (colour online) mobility edge trajectory vs doping. 
Circles (blue): band edges. Dash-dot line (red): Fermi en- 
ergy. Solid Lines (black): Mobility edges using actual DMFT 
conductivity. Dashed lines (red): Mobility edges using weak 
scattering approximation. Parameters: Ejr/t — —3,D/t — 6 

The mobility edge trajectories obtained by both meth- 
ods are plotted as a function of the doping in Fig [T] One 
finds that for either method, there are two mobility edges 
occuring symmetrically about the centre of the band, be- 
tween which lie the extended states. States lying between 
any one mobility edge and the band edge on the same side 
are localized, since it is the band tails which get localized 
first. For low doping values below x® MFT , all states in 
the band are localized, and both the mobility edges are 
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coincident at the band centre. As the doping increases, 
both the mobility edges proceed outwards from the band 
centre. Simultaneously, the bandwidth also increases 
(roughly as yfxD for low doping values), due to reduced 
scattering from polarons. The rate at which the mobil- 
ity edges move apart in energy with increasing doping is 
greater than the rate at which the bandwidth increases, 
so that ultimately, the mobility edges meet the corre- 
sponding band edges at a large enough doping (« 0.6 
for this parameter set). This means that beyond this 
doping value all states in the band are extended. This 
makes sense, because nearabout this large doping, there 
are no more ^-polarons left {W\ — 0). This means that 
the 6-electrons do not encounter any scattering, and the 
bandwidth is at its maximum. It is to be noticed that the 
mobility edge trajectories calculated by the two methods 
differ significantly only for small doping, while they co- 
incide for large doping. This is to be expected, since for 
smaller doping values, there are a substantial number of 
polaron scatterers, so that the weak-scattering assump- 
tion is not justified. 
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FIG. 2: Ejt vs x showing various phase boundaries, along 
with the DMFT parabola for comparison. 

We have also plotted the self-consistently determined 
value of chemical potential from the DMFT, on the same 
graph. It is found that even though the 6-band states 
begin to become occupied for x > x® MFT , but for a 
range of x thereafter, the occupied states are localized. 
The chemical potential proceeds towards the band centre, 
and meets the lower mobility edge at a doping x ^ nderson 
of the order of 0.4. Hence, for intermediate values of 

doping X® MFT < X < x Anderson^ S y S ^ em jg an A n _ 

derson Insulator, although DMFT predicts that it should 
be metallic. 

We now compare this result, obtained using effective 
field theories, with a more realistic, real space, finite-size 
simulation of the £ — b model (in the fully polarised fer- 
romagnetic phase at T=0), done by V. Shenoy et.al^. 
Details of this simulation will be found in Re^S. In this 
simulation, in addition to the Hamiltonian described be- 
fore, there is a long range Coulomb interaction between 
charges. This long-range Coulomb interaction (treated 



4 




i.2 0.25 0.3 0.35 0.4 0.45 0.5 

x (doping) 

FIG. 3: mobility edge trajectory vs doping from real space 
numerical simulation 

using Hartree-Fock approximation) prevents phase sepa- 
ration, and bring about a somewhat 'homogeneous' mix- 
ture of i and 6, as is assumed in DMFT. The presence of 
this long range Coulomb interaction, along with the com- 
pensating negative charge for the dopant atoms, gives an 
additional source of disorder that was not included in 
the effective medium calculations. However, this disor- 
der is correlated, as compared to the uncorrelated disor- 
der due to the £-polarons. One finds that this additional 
correlated disorder seems to make very little difference 
from DMFT results with regard to the 'critical' doping 
required for electrons to start occupying band states. 
This is observed in Fig [2j where the phase boundary 
for 6-state occupancy nearly coincides with the DMFT 
parabola given by Eg I15I 




1.2 0.225 0.25 0.275 0.3 0.325 
doping (x) 



FIG. 4: Doping variation of localization length evaluated at 
the Fermi energy (top) and at the center of the band (bottom) 

We refer to this 'critical' value of doping as x ° ccupancy . 
One finds that this x ° c ccupancy « xf MFT . However, in 
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this real space simulation, even after 6-band states be- 
gin to get occupied, the hole clumps within which b- 
electrons are allowed to hop (for on-site Udd — > oo) may 
not be percolating. We call the 'critical' doping for this 
percolation threshold as x P ercolatlon ; anc j this is higher 
than x ° ccu P anc v , Moreover, even within these percolating 
clumps, some of the 6-electron states may be Anderson 
localized. To check the character of these band states, 
we find the Inverse Participation Ratio (IPR), defined as 
follows: 



IPR 



(Ef hi 2 ) 2 



(16) 



where a* is the amplitude of the state at the i-th site. For 
an extended state, at goes as e lk ' Ri , so that \di\ — 1, and 
IPR — > -k — >• for large N. But for a state localized at 
the j-th site, a% — Sy, and IPR — > 1. Hence, IPR gives 
a good criterion for distinguishing extended states from 
localized ones. 

Since different realizations of the compensating dopant 
charge distribution give different set of states, one gets 
a distribution of IPR within each small energy range in- 
side the band. Hence, we discrctize the band into en- 
ergy intervals, and find the probability distribution for 
IPR within each such interval, by plotting histograms 
and corresponding frequency polygons, as discussed in 
detail in Ref^i. One finds from Fig [3] that the mo- 
bility edge trajectories once again proceed away from 
the band centre, and towards the edges as doping in- 
creases. Moreover, the chemical potential, determined 
from the band-filling in the same simulation, lies be- 
tween the lower band edge and the lower mobility edge 
for small doping. Hence, this simulation also predicts the 
existence of an Anderson Insulator phase for a certain re- 
gion of doping. Moreover, the value of x Anderson ~ 0.4 
obtained from Fig [3J also agrees closely with the effec- 
tive field theory calculations described before. Thus, for 
E JT = 3, the x fcu P ancy _ q.25, x P erc ° laU ° n w 0.33, while 

^.Anderson ^ q 

We have also calculated the localization length A as a 
function of doping using the PAW formalism. It turns out 
to be surprisingly small near x® MFT : namely of the or- 
der of a lattice spacing, even less for lower doping values. 



As one approaches x Anderson : \ a t Fermi energy diverges, 
as expected. A for states at the center of the band diverge 
for lesser x, since the mobility edge trajectories converge 
to the center at around x = 0.3. Variable range hopping 
is therefore expected to be observed in the Anderson In- 
sulator regime from x® MFT to x Anderson . The resistivity 
is therefore expected to go a s 23 ' 24 expiTo/T) 1 ^, where 
the temperature parameter Tq is given by Tq — 



k B N(E F ) ' 

a being the inverse localization length. From our calcu- 
lation of localization length, we estimate this parameter 
to lie between 10 7 — lO 8 ^ in this doping regime. 

In conclusion, we have extended the DMFT calcula- 
tions for the IMT within the two-fluid model 5 for man- 
ganites to incorporate localization effects. We find the 
existence of a prominent Anderson Insulator phase in ad- 
dition to the polaronic, percolative insulator phases seen 
before. Such an Anderson Insulator phase may well ex- 
plain the observations of variable range hopping (VRH) 9 
in manganites reported by experimentalists 2 ^—. The 
exceptionally small A in the Anderson Insulating regime 
should translate to a large value for the temperature pa- 
rameter for the VRH, as suggested by these studies. In 
the real space simulations by Shenoy et.al^, there exists 
a Coulomb glass phase with a soft Coulomb gap at low 
doping values prior to the occupancy of the 6-band states. 
This is due to the site-trapped i polarons, which also 
interact with each other by Coulomb interaction. This 
should give rise to a Shklovskii-Efros (SE-VRH) T 1 ' 2 
transport^!, which in this scenario will crossover into the 
Mott-VRH T 1 / 4 regime as the doping is increased, and 
the 6-band states begin to get occupied, before the fi- 
nal Insulator-Metal transition. Interestingly, both the 
SE-VRH and Mott-VRH regimes have been noticed and 
reported in the experimental literatur o 31 ^ 32 , and there ex- 
ists considerable debate between which one is more appli- 
cable to manganites. This is partly because of the techni- 
cal difficulty associated with distinguishing between var- 
ious power-laws in the exponent by fitting the resistivity 
data within a limited range of temperature 32 . Our cal- 
culation, together with the earlier ones by Shenoy et.ai, 
on the other hand, postulate the existence of both the 
SE-VRH and Mott-VRH regimes at different values of 
doping within the phase diagram for manganites. 
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